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ABSTRACT 



Various computational methods and operational computer codes used to predict and 
evaluate aerodynamic coefficients and flight performance of missile bodies are reviewed. 
Aerodynamic effects of symmetric and asymmetric flow separation are discussed, as are 
the differences inherent in estimating the properties of the resulting flowfields. The 
semi-empirical aeroprediction codes NSWC and MISSILE DATCOM are compared 
against experimental data for a variety of configuration geometries and flight conditions; 
the .MISSILE D.-\TCOM code is further used for a comparison with wind tunnel data 
for a Standard-type missile model. The NSWC and MISSILE DATCO.M codes are 
found to provide accurate prediction of normal force coefficients at both low and high 
angle of attack, although the nonlinear effects of separated flow are only partially cap- 
tured. Center of pressure coef['icients are generally underpredicted, but of the correct 
order of magnitude. The accuracy of drag coefficient prediction is seen to diminish as 
missile configuration geometr>’ becomes more complex. The NSWC program provides 
satisfactory prediction of pitch damping coefficients, while the MISSILE DATCOM 
output is inconclusive. The NSWC and MISSILE DATCOM aeroprediction codes arc 
considered suitable for preliminary design and aerodynamic analysis. 
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I. INTRODUCTION 



Flight stability and performance requirements are fundamental design consider- 
ations in tactical missile development. As missions and operational roles are redefined, 
so too are the aerodynamic conditions under which the missile must operate. A new 
flight environment may require only slight modifications to an existing system, or may 
involve significant design alteration and feasibility testing. The final result must, how- 
ever, continue to maximize the controllability and performance of the missile, thereby 
maintaining a dependable and effective weapons system. 

A. AERODYNAMIC MODELING 

The need for design alterations can arise from new employment methods, the pres- 
ence of new or upgraded threats, or an introduction of new technology. Recently the 
naval surface combatant community has seen the vertical launch system enter service 
aboard two new construction ship classes. Use of the vertical launch system for Anti-Air 
Warfare and Anti-Submarine Warfare places these weapons in an operational (light re- 
gime unforeseen during their development. Live-fire exercises of vertically launched 
weapons have been conducted, the results of which indicate that current design stability 
and performance require further investigation. During a test of a prototype thrust vector 
control (TVC) ASROC, the missile suddenly entered into an errant and unstable flight 
mode. Unable to control yaw plane motion, the weapon traveled a flight path nearly 90 
degrees off the intended down range trajectory’. Careful review of flight films and 
telemetry data indicates that an excessive aerodynamic yawing moment was induced 
through nose tip movement. Flight failure analysis suggests that this displacement in 
nose position most probably generated significant asymmetric vortex shedding, resulting 
in large out-of-plane moments, which continued even after symmetric reseating of the 
nose tip. [Ref 1). 

Although such asymmetric shedding of vortices can be induced by structural failure 
in the region about the nose, this effect is known to be a consistent flow phenomenon 
for missile bodies within the low Mach, high angle of attack regime. These conditions 
are typically encountered during the boost phase of a vertically launched missile. Many 
researchers have investigated the effects and properties of asymmetric flow as related to 
missile bodies. Experimental results indicate that flow separation commences at angles 
of attack of just a few degrees. At such relatively low incidence, the formation of lee side 
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vorticity occurs in a symmetric and largely steady manner. At angles of attack of 25 de- 
grees or more, asymmetric vortex development and shedding take place. While the 
asymmetric structure is nominally steady, a noted unsteady behavior is observed to oc- 
cur at angles of attack near 55 degrees [Ref 2]. Under certain flow conditions, the vortex 
pattern is known to rapidly fluctuate between nearly symmetric and highly asymmetric 
structures [Ref 3j. The exact mechanism of asymmetric vortex generation about missile 
body configurations is not yet fully understood, although the general consensus of in- 
vestigators links this phenomenon to Reynolds and Mach numbers, nose fineness ratio, 
and the angle of attack [Ref 4]. It is interesting to note that asymmetric modeling usu- 
ally requires the introduction of a perturbation about the nose region, thus creating 
conditions very similar to those experienced by the TVC ASROC. 

The rapid application of new technologies towards improved capability and newly 
configured missiles demands a responsive method for the evaluation of design modifica- 
tion proposals. Such performance evaluations must include a comprehensive analysis of 
the dynamic factors which exist throughout the operating envelope of the weapon: in- 
creased maneuverability, high structural loads, high angle-of-attack aerodynamic effects, 
control surface limitations, and Mach number effects. By accurately modeling these 
conditions, the forces on and the moments about the missile body can be predicted and 
the flight trajectory estimated. A proper and complete description of the missile-flowfield 
interaction is both necessary and complex, but often must be approximated in a manner 
which permits a simplified application while still fulfilling the scope of the research ef- 
fort. Such aeropredictive models rely on iterative mathematic computations and data 
base comparison, and are well suited for computer system implementation as a coded 
program. 

A computer-driven solution can be based entirely on the theoretical relationships 
and derived equations that describe a body in motion through a fluid. Codes con- 
structed in this manner are extremely complex, but may offer a precise numerical pro- 
cedure for solution. Such numerical programs require a great amount of computer time, 
and remain largely within a research stage. Many of these numerical programs employ 
new modeling techniques or innovative computational schemes. Less complicated codes 
combine computations and relevant empirical results. There is a wide range in both the 
capability and modeling approach of codes within this category. These programs gener- 
ally show a good degree of consistency with accepted experimental data or alternate 
codes. 
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The accuracy of any code is, at best, limited to the level with which the missile- 
flowfield interface has been modeled. A purely mathematical approach is conditionally 
capable of a complete solution, but may require exorbitant development costs and re- 
strictive computer time. Furthermore, as these theoretical codes are somewhat tailored 
to specific applications, they may be sensitive to initial conditions and initial data inputs. 
These characteristics may make such a theoretical code untenable and undesirable for 
general purpose use. Semi-empirical methods seek a valid and accurate combination of 
computational mechanics and experimental observation, and are frequently developed 
on and extended to the smaller computer systems normally associated with research and 
academic institutions. Many semi-empirical programs currently in use are accurate to a 
level satisfactory for preliminary or intermediate design. While these codes are less ca- 
pable than theoretical programs, they are easier and faster to operate, and afford the 
user a wide range of input options. Semi-empirical codes pennit an accelerated prelimi- 
nar}’ design process, and can greatly reduce the amount of costly wind tunnel testing 
required during this phase. In view of the flexibility and efficiency to be gained, 
aeroprediction modeling codes have become essential tools for aerodynamic research and 
development. 

B. OBJECTIVES 

The objectives of this thesis are to research the available operational predictive 
missile codes and identify those suitable for installation and operation on the mainframe 
computer system at the Naval Postgraduate School, in support of the Weapons Engi- 
neering academic curriculum and current high angle of attack missile research. Acquired 
programs are to possess capabilities for modeling conventional missile configurations in 
both subsonic and supersonic Mach regimes. Force and moment prediction is to include 
analysis at high angle of attack, with the desired inclusion of asymmetric separation 
features. Predicted quantities for selected input data are to be compared and discussed 
for each code. Finally, simulations based on the geometrx' of a Standard-type missile 
model will be made and compared against Naval Postgraduate School wind tunnel data 
of this model, the testing of which supports an ongoing research effort for the Naval 
Surface Warfare Center at Dahlgren, Virginia. 
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II. PREDICTION PROGRAMS FOR AERODYNAMIC ANALYSIS 



The capability and validity of a computational prediction code are primarily deter- 
mined by the detail and depth to which the missile-flowfield interactioYi has been mod- 
eled. A program typically becomes more involved as consideration is given to complex 
and non-traditional configuration geometries, and is further compounded for application 
of nonlinear aerodynamic effects. The development of a satisfactory’ model cannot, 
however, be done irrespective of computer system requirements. The construction of an 
aeroprediction code should include the following criteria: 

• flexibility as to allowable configuration geometries. 

• applicability to a sufficiently broad range of flight conditions. 

• accuracy at a level congruent with development objectives, research goals and 
operational requirements. 

• run time and computer s> stem requirements reasonably available to the intended 
user base. 

Neglect in consideration of these factors may severely limit the operability of a program 
or render it altogether unsatisfactory. 

Most aerodynamic prediction codes can be categorized as either theoretical or 
semi-empirical, and are discussed in this manner. A brief review of several techniques 
and modeling principles illustrates the fundamental dilTerences of various methods, as 
well as similarities inherent to missile prediction codes. Selected operational codes are 
presented for comparison along these lines. A review of missile aerodynamics and oper- 
ational prediction programs is given by Lacau [Ref 5). 

A. THEORETICAL PREDICTION CODES 

Codes of this type perform a purely computational analysis based on numerical 
methods. Such programs model the governing partial differential equations (based on 
certain simplifying assumptions) by finite mathematics. These codes have shown excel- 
lent results when applied to three dimensional missile aerodynamics, and have provided 
investigators with valuable information pertaining to local flowfield characteristics and 
complex flow mechanisms [Ref 6]. Although severe computer constraints currently limit 
the application of such programs, rapid advances in computer technology make the 
computational fluid dynamics code a promising approach for future work. Numerical 
methods of high incidence missile aeroprediction are developed using the Navier-Stokes 
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equations, the Euler equations or the linearized potential flow equation (Prandtl- 
Glauert). 

1. Navier-Stokes Methods 

The Navier-Stokes equations should be capable of describing any flowfleld over 
a missile body under any conditions. Predicted quantities include the effects of separated 
flow' and rapidly fluctuating asymmetric vorticity. With veiy few exceptions, simplifi- 
cations must be made to the full Navier-Stokes equations due to limitations in computer 
system resolution and an incomplete understanding of the physics of turbulence. 
Normally, the fluctuating components are time-averaged to yield the Reynolds-averaged 
equations. The Reynolds equations can be extended to the most complex flow condi- 
tions, however, a suitable turbulent flow' model must be incorporated to provide closure 
of the solution, such as the algebraic method of Baldwin and Lomax [Ref 7]. Ihe 
Reynolds averaged equations can be further reduced to the thin-layer (neglecting 
streamwise viscous terms) or parabolized (neglecting unsteady terms and streamwise 
viscous diffusion) Navier-Stokes form. These equation types can be solved using time- 
marching and space-marching techniques for either turbulent or laininar flow' conditions. 
A selection of Navier-Stokes codes is presented in Table 1 on page 6. 
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Table 1. NAVIER-STOKES MISSILE CODES 



CODE NAME 


CODE TYPE 


LAMINAR 

TURBULENT 


MACH REGION 


GEOMETRY 


ARC3D 


FULL 


LAMINAR 

AND 

TURBULENT 


SUBSONIC 

TRANSONIC 

SUPERSONIC 


CONVENTIONAL 
OR ARBITRARY 


F3D 


THIN LAYER 


LAMINAR 

AND 

TURBULENT 


SUBSONIC 

TRANSONIC 

SUPERSONIC 


CONVENTIONAL 
OR ARBITRARY 


NASBMG 


FULL 


LAMINAR 

AND 

TURBULENT 


SUBSONIC 

TRANSONIC 

SUPERSONIC 


CONVENTIO.NAL 
OR ARBITRARY 


LWIN 


FULL 


LAMINAR 


SUBSONIC 

TRANSONIC 

SUPERSONIC 


CONVENTIONAL 
OR ARBITRARY 


HNS 


PARABOLIZED 


LAMINAR 


MACH > 1 


CONVENTIONAL 
OR ARBITRARY 


PNSFVM 


PARABOLIZED 


LAMINAR 


MACH > 1 


CONVENTIONAL 
OR ARBITRARY 



Source: [Ref. 5, p. 1-43] 



or the relatively few Navier-Stokes codes in operation, one frequently applied is F3D, 
developed by Steger, Ying and SchilT[Rer S]. A thin layer, time accurate code, F3D has 
been used by Degani and Schiff to study three dimensional subsonic flow about a slender 
body of revolution at high angle of attack [Ref. 9]. For laminar flow, the coefficient of 
viscosity is taken from Sutherland's law, while turbulent conditions make use of an 
eddy-viscosity model introduced by Degani and SchifF[Refs. 10,11]. Further work has 
been done regarding vortex unsteadiness and the effects of spatial disturbances on vortex 
asymmetry at large incidence [Refs. 12,13]. Degani and Schiff report generally satisfac- 
tory computational results from F3D compared to experimental data, with computation 
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limes on the order of 30 seconds per iteration on a Cray supercomputer. Similar work 
has been done using an incompressible thin-layer Navier- Stokes code, FMCl. Hartwich 
and Hall, using an eddy-viscosity model based on the two-layer, zero-equation form of 
Baldwin and Lomax, have applied the FMCl program to determine leeward pressure 
distributions on a tangent ogive body with large crossflow separation [Ref 14]. Numer- 
ical solutions are good to excellent in comparison to experimental data, and represent 
an improvement to the results obtained using the F3D code with the turbulent flow 
model of Degani and SchifT. The required CPU times (per grid point per iteration at a 
Reynolds number of 2.0E5) are roughly 29 microseconds on a CDC Cyber 205 machine 
and 58 microseconds on a Cray 2. For low Mach regime analysis, the FMCl code offers 
greater efficiency and increased accuracy over the F3D program, which appears to be 
characteristic of the incompressible flow codes. [Refs. 15,16.17]. 

2. Euler Methods 

Euler equations represent the Navier-Siokes equations in which the viscous and 
conduction terms are ignored. Euler equations are descriptive of inviscid rotational or 
irrotational flow independent of Mach number, and can be applied to flow conditions 
with shock formation and vortex sheets. The steady equations are hyberbolic type partial 
dilferential equations, and are solved through space-marching techniques which restrict 
the application to supersonic speeds. The procedure for the solution of the unsteady 
equations is to advance the complete flow variable array in time until convergence oc- 
curs at some asymptotic limit; this approach is valid for both subsonic and supersonic 
Mach regimes, but neglects fluctuation terms induced by movement of the missile and 
associated flowfield changes. Various operational Euler codes are presented in Table 2 
on page 8. 
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Table 2. EULER BASED MISSILE CODES 



CODE 

NA.ME 


EQUATIONS 


METHOD 


NUMERICAL SCHEME 


MACH 




STEADY 

OR 

UNSTEADY 


CONSERVATIVE 
OR NON- 
CONSERVATIVE 


FINITE VOLUME 
DIFFERENCES 
OR ELEMENTS 


IMPLICIT 

OR 

EXPLICIT 


CENTERED 
OR UN- 
CENTERED 


< 1.0 


>1.0 


EUFLEX 


U 


C 


FV 


BOTH 


UC 


YES 


YES 


ELLBMG 


U 


C 


FV 


E 


C 


YES 


YES 


EULER3D 


U 


C 


FV 


E 


C 


YES 


YES 


EL'LSSM 


s 


C 


FV 


I 


C 


NO 


YES 


FLL3C 


u 


C 


FV 


E 


UC 


YES 


YES 


MISSILE 


s 


C 


FD 


E 


C 


NO 


YES 


MLSE 


s 


C 


FD 


E 


C 


NO 


YES 


SANDI.VC 


s 


BOTH 


FD 


E 


UC 


NO 


YES 


SWINT 


s 


C 


FD 


E 


C 


NO 


YES 


WING2A 


u 


c 


FV 


E 


C 


YES 


YES 


ZEES 


s 


c 


FD 


E 


UC 


NO 


YES 



Source: [Ref. 5, p. 1-44] 



While not applicable to viscous flow analysis, Euler codes are powerful 
aeroprediction programs. Currently operational codes are capable of describing the flow 
conditions about any configurational geometry, including separated flow and shock for- 
mation effects. Two such codes are SWINT and ZEUS. These programs provide excel- 
lent accuracy through the use of finite difference solutions. Operator use is simplified in 
that a minimal amount of preliminary set up is required. Described as robust, both 
SWINT and ZEUS are receiving extensive use in missile aerodynamic prediction and 
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design proposal performance evaluations. Techniques used in the structure of these 
codes, primarily the Godunov method of solution of the Euler equations, show partic- 
ularly good results. Current work is in progress to expand such techniques to include 
viscous effects through application to Navier-Stokes equations. [Refs. 18,19,20]. 

Although reduced in complexity from the Navier-Stokes form, Euler codes re- 
quire computer times on the order of comparable Navier-Stokes programs, and are 
normally not a viable method of conducting general aerodynamic analyses. A particular 
limitation exists in the solution of separated flow along a smooth surface, which requires 
the determination of cross flow separation lines for application of the Kutta condition. 
As appropriate nonlinear data bases are frequently unavailable, this information must 
be obtained through boundary layer flow solutions or from a suitable Navier-Stokes 
code. 

3. Potential Flow Methods 

The linearized potential flow equation provides a method to describe flow con- 
ditions with induced perturbation velocities slightly different than those in the 
freestream. Codes of this type are no longer purely numerical, reflecting instead an ap- 
proximate mathematical form with certain simplifying assumptions required by the 
modeling approach. Most aeroprediction applications of the linearized potential 
equation are centered about a panel method. Panel methods are general numerical 
techniques used to describe Oowfields about arbitrary bodies. While a large number of 
different applications have been developed, the vortex panel method is frequently used 
to evaluate the aerodynamic forces induced through missile-flowfield interaction. 

Consideration is given to an arbitrary three dimensional body within an inviscid, 
incompressible How. Vortex filaments are modeled as a vortex sheet, and the body sur- 
face is conceptually wrapped within this sheet. The approach is to explicitly solve for the 
vortex sheet strength (per unit length) such that the continuous conform of the body 
surface satisfies both the streamline conditions for the surrounding flow and the Kutta 
condition at the points of crossflow separation. In the absence of a closed form analytic 
solution, a numerical solution is obtained through the reduction of simultaneous 
equations. This requires approximating the vortex sheet as a series of panels about the 
body. The unknown vortex strength is taken as constant across each panel in first order 
solutions; higher order techniques model the panel strength functionally, such as a linear 
variation with surface displacement. The midpoint of each panel is selected as a control 
point at which boundar}' conditions of zero velocity of the crossflow normal component 
are enforced. The Prandtl-Glauert equations are used to determine the panel induced 
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velocity at any point within the flow. Taken about each control point, a geometn.' de- 
pendent set of simultaneous equations is obtained which relates crossflow velocities to 
local panel vorticity. These equations, coupled with both the Kutta and streamline con- 
ditions. are solved to yield individual vortex sheet strengths. Flow tangency conditions 
are applied to determine the local variations in the flowfield above tli6 sheet, which al- 
lows computation of pressure distributions and induced forces on the body. 

The above procedure presents a valid scheme for the prediction of aerodynamic 
quantities on an arbitrary body in subsonic flow. A severe limitation exists, however, in 
that the linear nature of the solution restricts analysis to inviscid flow at low angles of 
attack. As such, application to separated flowfields is precluded due to the nonlinear 
mechanics of flow separation and vortex formation. The occurrence of flow separation, 
particularly at high angles of attack, causes a dramatic variation in both the strength and 
distribution of vortices, which in turn produces a wide fluctuation in the induced aero- 
dynamic forces and moments. In view of the previously referenced firing exercises, the 
analysis of such flight conditions is obviously of considerable interest. 

\\'ork has been conducted in an attempt to extend panel methods to include 
nonlinear eflects arising from vorticity and compressibility [Ref 21]. Paneling techniques 
used to modify the linearized potential codes PFPl through PFP5 have been applied by 
Van Tuyl for calculations about a missile at high angle of attack. The vortex wake is 
appro.ximated by vortex sheets attached along pre-determined separation lines, with flow 
tangency conditions satisfied by the inclusion of source panels. Leeward pressure dis- 
tributions are found to be in good agreement with thin-layer Navier-Stokes calculations 
and experimental data. The windward prediction results arc inconclusive, with the values 
in good agreement with the empirical data base, but not with the Navier-Stokes sol- 
utions. [Refs. 22. 23 ]. Mendenhall and Perkins have conducted significant studies into 
the vortex-induced characteristics of missiles within nonlinear and unsteady flow condi- 
tions. This work, using vortex lattice panel methods and vortex cloud modeling, has 
shown that nonlinear eflects can be captured using panel techniques, although the re- 
sulting code is quite complex and expensive to develop and run. The prediction code 
NOZVTX was found to accurately describe the separated flowfield to the Iccside of a 
missile at a level suitable for engineering or preliminary design. It is interesting to note, 
however, that predicted values of windward pressure distribution show close agreement 
with experimental data, but are lacking to some degree in comparison with Navier- 
Stokes solutions, as was similarly reported by Van Tuyl. [Refs. 24,25]. 
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The development of aeroprediction codes based on the linearized potential flow 
equations represents a departure from the purely theoretical structure of Xavier-Stokes 
and Euler methods. In constructing these codes it becomes necessary to make critical 
assumptions and modeling approximations, such as in the distribution of panels and 
vortex strengths, the configuration of body contour relative to the freestream, and the 
determination of crossflow separation points. Errors in these input criteria may well re- 
sult in distorted output to a point that meaningful interpretation is inhibited. While these 
codes are appreciably less involved than those derived from the wholly mathematical 
description, their approximate theoretical form is less precise. The results of these pro- 
grams have been excellent, however, and certainly warrant further development. Al- 
though currently considered an advanced design tool, these codes will become 
increasingly useful as computer system capabilities continue to advance. A selection of 
linear potential equation codes is shown in Table 3 on page 12. 
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Table 3. LINEARIZED POTENTIAL BASED MISSILl 


E COD 


ES 


CODE NAME 


VORTEX SEPARATION 


MACH 


GEOMETRY 


DEMON/LRCDM 


BODY AND FIN VORTEX SEPAR.ATION 
VORTEX MODEL WITH PANEL METHOD 


M > 1 


CONVENTIONAL 
OR ARBITRARY 


DM31.nl 


FIN VORTEX SEPARATION WITH 
TRAILING EDGE WAKE RELAXATION 
VORTEX MODEL WITH PANEL METHOD 


M > I 


CONVENTIONAL 
OR ARBITRARY 


HOP 




M < 1 
M > 1 


CONVENTIONAL 
OR ARBITRARY 


\A.\C 




M > I 


CONVE.NTIONAL 
OR ARBITRARY 


NFKRSUB/SUP 


BODY AND FIN VORTEX SEPARATION 
VORTEX MODEL WITH PANEL METHOD 


M < 1 
M > 1 


CONVENTIONAL 


NLRAERO 




.\1 < I 
M > 1 


CONVE.NTIONAL 
OR ARBITRARY 


NWCDM/NSTRN 


BODY AND FIN VORTEX SEPARATION 
WHH TRAILING EDGE WAKE RELAXATION 
VORTEX MODEL WITH PANEL METHOD 


M > 1 


CONVENTIONAL 


PANAIR 




M < 1 
.\I > 1 


CONV E.NTIONAL 
OR ARBITRARY 


PANEL 




M < 1 


CONVENTIONAL 
OR ARBITRARY 


QRFL/DEMONI 


BODY AND FIN VORTEX SEPARATION 
VORTEX .MODEL W ITH PANEL METHOD 


\I < 1 
M > 1 


CONVENTIO.NAL 


MBC 


BODY AND F IN VORTEX SEPAR.ATION 
VORTEX MODEL WITH PANEL METHOD 


M < 1 


CONVENTIONAL 
OR ARBITRARY 



Source; [Ref. 5, p. 1-45] 



Theoretical programs olTer the flexibility of unrestricted geometry configura- 
tions and the fundamental ability to accurately describe even the most complex missile- 
flowfield interactions. In addition to computational aerodynamic predictions, the 
modeling techniques ofTer important explanatory' insight into the mechanics which exist 
for nonlinear flow. The employment of such codes is currently restricted to research fa- 
cilities which possess the tremendous computer systems required, and is generally not a 
feasible option for basic design work or aerodynamic prediction. Until computer tech- 
nology delivers these codes into wide distribution and commonplace use, smaller re- 
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search institutions and academic concerns must rely on alternate means to conduct 
efficient and rapid preliminary’ evaluations of stability and performance. 

B. SEMI-EMPIRICAL PREDICTION CODES 

Semi-empirical codes are primarily designed to provide for rapid and efficient aero- 
dynamic analysis of traditional or conventional configurations. These programs attempt 
to combine approximate theoretical relationships and empirical data such that the out- 
put agrees with observed experimental results. By excluding routines which attempt to 
interpret various flow parameters, such as the nonlinear mechanics of the high angle of 
attack regime, these methods can be greatly reduced in complexity as compared to the- 
oretical models. Semi-empirical codes are intended for practical, general purpose anal- 
ysis, and can normally be operated on computer systems considerably less capable than 
those required for numerical codes. The principal methods of semi-empirical codes are 
slender body and linearized potential flow theory, which are most often utilized in a 
component build-up approach. 

The component build-up procedure calculates the individual aerodynamics of each 
elementary portion of the missile body, estimates the effects of component interaction, 
and attempts to synthesize a solution based on the total body. While this method is 
generally restricted to conventional configurations, a component build-up approach can 
offer significant reductions in development and usage costs, good to excellent accuracy 
for traditional missile geometries, easy extension to parametric analysis, and simple 
inclusion of experimental data [Ref 26]. As the alternative to component build-up is a 
panel method geometry description, the majority of preliminary' design missile codes are 
structured around a build-up procedure. 

Semi-empirical aeroprediction programs can generally be separated into two distinct 
categories. One approach is an empirically oriented solution based on correlation to 
functional data for various component geometries, flight conditions, and angles-of- 
attack. The remaining technique is a more involved computational method based on 
crossflow modeling with the addition of body vortex theory. Several of the fundamental 
methods used are discussed below; a more complete background of empirical and semi- 
empirical routines has been prepared by Wardlaw [Ref 27]. Various empirical and 
semi-empirical codes are presented in Table 4 on page 14. 
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Table 4. EMPIRICAL AND SEMI-EMPIRICAL MISSILE CODES 



CODE 

-NAME 


-MACH 

RANGE 


INCIDENCE 

RANGE 

(DEGREES) 


ROLL 

ANGLE 

(DEGREES) 


AXIAL 

FORCE 


STATIC 

STABILITY 


DYNAMIC 

derivatives 


ABACUS 


M < 5.0 


-90 to 90 


-180 to 180 




CN.Cm.Cn 

Cy.Cl 


Cmq.Cmad 

Cnr.Clp 


AERA-M 


I.2<-\I<5.0 


0 to 30 


0 10 45 


CA.CAw 

CAf.Cab 


CN,Cm 


Cmq.Cmad 


BAKER 


0.6 <M <3.0 


0 to 180 


0 




CN,Cm 




CASAERO 


-M <4.0 


0 to 30 


0 




CN.Cm 


Cmq.Cmad 

Cnr.Clp 


DORRAM 


M <4.0 


0 to 90 


0 to 360 


CA.CAw 

CAf.Cab 


CN,Cm 

Cy,Cn 


Cmq.Cmad 

Cnr.Clp 


MAP 


M<8.0 


0 to 15 


0 to 90 


CA.CAw 

CAf.Cab 


C.\,Cm 

Cy,Cn,Cl 




MISSILE 

DAICOM 


M<8.0 


0 to 90 


0 to 360 


CA.CAw 

CAf.Cab 


C\,Cm 

C>,Cn,Cl 


Cmq.Cmad 

Cnr.Clp 


MISSll E 
ON ERA 


-M <4.0 


0 to 25 


0 to 360 




CN,Cm 

Cy,Cn,Cl 




MISSILE I 


M<4.0 


0 to 45 


0 to 360 


CAvv.CAf 


CN.Cm 

Cy.Cn.Cl 




MISSILE 2 


M < 5.0 


0 to 45 


0 to 90 




CN.Cm 

Cy.Cn.Cl 




MISSILE 3 


M <4.5 


0 to 45 


0 to 90 




CN.Cm 

Cy.Cn.Cl 




-NSWC 


M <8.0 


0 to 45 


0 to 360 


CA.CAw 

CAf.Cab 


CN.Cm 

Cn.Cl 


Cmq.Cmad 

Cnr.Clp 


S/IIABP 


M>2.0 


-90 to 90 


0 to 360 


CAw.CAf 


CN.Cm 

Cy.Cn.Cl 





Source: [Ref. 5, p. 1-42J 



1. Empirical Methods 

Empirical techniques for aerodynamic analysis are often simplified models of 
some frequently observed and well documented conditions. These routines attempt to 
replicate the experimental results through the proper use of data base information. 
Perhaps the most well known method for body alone aerodynamic prediction is the so 
called crossflow analogy of Allen and Perkins [Refs. 28,29). In this method, the flow over 
a body at incidence is treated as the resultant of normal flow and crossflow components. 
The body normal force represents the inviscid, or potential flow contribution, and is 
calculated using slender body theory. The viscous crossflow term is evaluated through 
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analogy to the flow about a cylinder, in which the crossflow drag coeflicient of the body 
is equated to the steady state drag on the cylinder. Further work by Jorgensen has ex- 
tended this method to transonic regimes for angles of attack up to 180 degrees. 
[Refs. 30,31]. Another empirical approach is the impulsive flow analogy of Kelly [Ref 
32]. In this method the crossflow about an inclined body is likened to the flow around 
an impulsively started cylinder. Specifically, the assumption is made that the crossflow 
plane is swept uniformly down the length of the body. The translation of the crossflow 
plane is viewed analogously to the leeside flowfield produced behind the inclined cylin- 
der. Flowfield properties induced by the changing nose radius of the body are captured 
by allowing the cylinder to expand [Ref 27, pp. 3-8]. Further work using this method 
has been done by Thomson in which an updated data base for impulsive cylinders 
(Sarpkaya) was available [Refs. 33,34]. Many variations of these models have been de- 
veloped which incorporate adjustments so as to include the effects of more complex 
phenomena such as those related to boundar>’ layer flow, base regions and nose 
bluntness. 

A widely used operational code is MISSILE DATCOM of McDonnell Douglas 
[Refs. 26,35]. This program uses component build-up techniques to provide aerodyna- 
mic prediction for traditional missile configurations and arbitrarx’ geometries. Static co- 
eflicient and dynamic derivative computations can be had for bodies at high angles of 
attack, but no analysis of asymmetric flow separation is possible. For conventional 
symmetric missile bodies, the effects of angle of attack are estimated using an extension 
of the .Allen and Perkins method, as expanded by Jorgensen. Potential flow solutions are 
taken from the second order shock equations (SOSE), and Van Dyke hybrid theory; the 
viscous crossflow solutions are taken from correlation with empirical crossflow drag data 
of -Messersmitt-Bolkow-Blohm. Configuration synthesis includes vortex tracking and 
strength calculations for body and lifting surface vortices. The development of the 
.MISSILE DATCO.M source code was largely driven by an extensive review and com- 
parison of existing data bases and computational methods. This program has shown 
simple operation on a number of different computer systems, with compilation times on 
the order of several CPU seconds for CDC Cyber machines. While the accuracy of spe- 
cific quantities is dependent on input geometry and flight conditions, the overall accu- 
racy of the code is deemed to be at a level suitable to support general purpose, 
preliminary design studies. 

,A similar aeroprediction code has been developed by Devan, w’hich will be re- 
ferred to here as NSWC [Ref 36]. Work on this program continues as an ongoing re- 
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search project of the Naval Surface Warfare Center (NSWC), with the source code 
presently in a fourth stage of development. Current efforts are directed at extending the 
range of application to higher Mach and higher angle of attack flight regimes. Prediction 
accuracy for dynamic derivatives is being improved as well. The NSWC aeroprediction 
code is structured as a component build-up approach for axisymmetric, traditional mis- 
sile and projectile geometries. This code offers the user a wide variety of options, but 
requires a more detailed input geometiy’ than does MISSILE DATCOM. High angle 
of attack aerodynamics do not include consideration of asymmetrically separated flow, 
and are restricted to flight conditions above the subsonic regime. Crossflow modeling 
follows the Allen and Perkins approach. Inviscid flow solutions are taken from SOSE, 
Tsien first order methods, Wu and Aoyoma computational methods, and solution data 
for Euler equations. High angle of attack viscous solutions are derived from empirical 
crossflow data which are maintained on a separate tape or input file. Aerodynamic 
quantities are predicted with the accuracy required for preliminar}' or intermediate design 
analysis and performance evaluations. Compilation times (per reference flight condition) 
are in the CPU seconds for a CDC 865 Cyber machine. 

2. Senii-Einpirical Methods 

Semi-empirical techniques attempt to model the flow conditions within the 
crossflow plane. .Many of these models use potential flow theor\', based on the impulsive 
flow description of leeside separation and trailing wake. As previously noted, the non- 
linear potential equation can be reduced through the application of perturbation princi- 
ples. If a slender body approach is assumed, the flow solution in the crossflow plane may 
be taken as incompressible which further reduces the linear velocity potential equation 
to a Laplacian form. The resulting linear partial dilTerential equation is significantly 
easier to solve than the nonlinear full potential equation, but such linear, inviscid sol- 
utions are no longer exact and are not suitable for analysis in nonlinear regimes, i.e., 
transonic flight speeds at high angles of attack. Semi-empirical methods generally com- 
bine computational routines and empirical data with which the effects of leeside sepa- 
ration and vorticity can be accurately modeled. 

Early semi-empirical techniques to analyze separated, nonlinear flow, made use 
of an impulsive flow model developed by Bryson [Ref 37]. In this approach, point 
vortices are superimposed on the potential flow solution about a cylinder. Vortex wakes 
are assumed to roll up into concentrated filaments, which remain attached to the body 
by feeding sheets. Vortex strengths are determined by application of the Kutta condition 
along empirically determined separation points. The motion of vortices along the body 
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length is such that the net force acting on the sheet and filament is zero, simulating an 
equilibrium condition between induced lift of the vortex and the pressure distribution 
across the sheet. Asymmetric separation can be modeled by introducing an initial per- 
turbation to the vortex positions about the nose, although the manner in which this 
perturbation is applied is apparently a strong determinant in the properties of the re- 
sulting asymmetric fiowfield. This approach has been further extended to body vortex 
models, such as in the work of Schindel and Wardlaw. [Refs. 38,39]. 

Other crossflow modeling work on viscous, unsteady flows has been conducted 
using multivortex or vortex cloud techniques. Numerous point vortices are used to rep- 
resent each wake vortex in the leeside crossflow plane. The developing vortices are peri- 
odically introduced into the fiowfield in the vicinity of crossflow separation points, and 
are modeled to roll up into large clusters of concentrated vorticity. Assuming that each 
point vortex maintains a constant strength after formation, vorticity transport equations 
can be taken from the Milne-Thomson circle theorem which allows calculation of vortex 
movement within the leeside fiowfield [Ref 40]. This approach was used by Angelluci to 
study symmetrically separated flows, and extended by Wardlaw to include asymmetric 
conditions [Refs. 41, 42 ]. Perhaps the most complete work is the discrete vortex cloud 
model of Mendenhall [Ref 43]. Improved methods for vortex tracking and strength 
calculations have been developed by Nielsen, et al., which are frequently adapted for use 
in currently operational codes [Refs. 44,45]. 

An early semi-empirical method was developed by Tidier and Bateman 
[Ref 46]. The basic approach couples a vorticity-conservation technique with a poten- 
tial crossflow model. The effects of incidence on vortex strength and location are based 
on experimental data and theoretical results for both symmetric and asymmetric condi- 
tions. The model considers the strength and spacing of growing and shedding vortices. 
The first two vortices to form behind the nose are treated individually, with strength 
contributions from potential flow about the nose. As indicated by flow surveys, the sol- 
ution method assumes that the remaining vortices are nearly equal in strength in relation 
to the boundary layer separation about the body. Solution of the first two vortices is 
taken from potential flow theory, while the so-called "street" vortices are evaluated using 
data on vortex strength as a function of angle of attack. This code achieves a fair degree 
of accuracy in the prediction of static coefficients and dynamic properties, although the 
onset of asymmetric vorticity is less satisfactory. Nielsen and Smith have developed a 
comprehensive and accurate model using similar methodology to describe vortex wake 
distributions for bodies at incidence [Ref 47]. 
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Codes which contain such semi-empirical computational schemes are normally 
developed for use in conjunction with panel method techniques. The mathematical 
routines require significantly increased computer memor>' and run time in comparison 
to more empirically oriented codes. Thus, the majority of such body vortex codes re- 
main unsuitable for general research or preliminary design studies. An exception to tliis 
loose rule may be the VORSTAB 11 aeroprediction code. This semi-empirical code is 
primarily intended for fighter aircraft configurations, but includes the capability for 
missile geometry analysis. Both symmetric and asymmetric forebody separation are cal- 
culated using slender body theory' and free vortex filament modeling. While symmetric 
vortex separation has been previously evaluated using vortex lattice techniques, the 
ability to accurately model asymmetric separation as an additional boundary' value sol- 
ution to a slender body, discrete vortex problem has only recently been demonstrated 
by Lan and Chin [Ref 48]. To avoid the necessity of turbulent model Navier-Stokes 
solutions, boundary layer separation effects are evaluated through the use of nonlinear 
section data as introduced by Lan [Ref 49 |. This is accomplished in a iterative match- 
ing procedure of nonlinear section lift data and lifting surface theory. Although not 
currently in distribution, discussions with the investigators indicate that this very robust 
code is suitable for use on most computer systems, with approximately two megabytes 
(.MB) required for an average compilation. [Ref 50]. The dramatic reduction in com- 
puter system requirements to operate this highly capable code should significantly in- 
crease the level of aerodynamic research and design which can be routinely undertaken 
by both academic and industrial concerns. 
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III. SELECTION AND INSTALLATION 



The selection process for aeroprediction codes was primarily influenced by the fol- 
lowing factors: availability of the source code and documentation, compatibility of 

source code structure, language and operation time requirements generally consistent 
with the computer center facilities of this school, and a demonstrated capability for 
treatment of traditional missile configurations at high angles of attack. It was addi- 
tionally desired that at least one code possess the ability to model asymmetric vortex 
shedding. 

The availability of prospective programs was examined as regards to ownership of 
source coding, that is, public domain information versus reserved proprietary' rights. 
Operational missile codes developed under federal (Defense Department) funding are not 
always public domain: innovative routines, data sets or modeling techniques may remain 
the property of the facility actually preparing and researching the source code. Proprie- 
tary source coding, if available at all. may have to be purchased, and may include re- 
strictions on application and release. 

The mainframe computer system at the Naval Postgraduate School is an IBM 
3033 4381 Network. Operations exist to support high volume batch processing and 
general terminal timesharing for the students, faculty, stafT and tenet commands. Ter- 
minal interactive use is performed by the two CPU 3033AP system, with 16 MB 
processor storage. Batch processing is accomplished by the 3033U and 43S1-M1 combi- 
nation. Programming languages include VS Fortran and Watfor 77. In addition to the 
IBM mainframe. VAX and micro- VAX workstations are available for use. 

A. SOURCE CODE PROCUREMENT 

In reviewing the literature concerning operational codes and computational meth- 
ods, it became apparent that numerical programs and paneling techniques were not vi- 
able options for use with the available computer system. Nor were such codes 
particularly desired, in that the intended application is largely in support of introductory 
level aerodynamic research. As the majority of semi-empirical missile codes are roughly 
equivalent in capability, selection was more heavily influenced by availability. 

In addition to the NSWC and .MISSILE DATCO.M codes, consideration was given 
to the data base oriented .MISSILE 1, .MISSILE 2 and .MISSILE 3 codes of Nielsen 
Engineering and Research [Refs. 51, 52 ]. In that the .MISSILE (1,2,3) codes provide 
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little additional aerodynamic prediction, excepting an increase in allowable control de- 
flection, the NSWC and MISSILE DATCOM programs were selected. The KSWC 
source code and documentation were obtained from Mr. Michael Armistead (G23), 
Naval Surface Warfare Center, Dahlgren, Va. The MISSILE DATCOM source code 
was provided by Mr. William Blake, Flight Dynamics Laboratory (AFWAL/FIGC), 
Wright-Patterson Air Force Base, Oh. Both of these codes have been used extensively 
to study the aerodynamics of missile bodies, and can be operated on the IB.M 3033 sys- 
tem. Neither of these two codes, however, is capable of performing analysis of asym- 
metric vortex shedding at high angles of attack. 

Identification of a suitable, asymmetric capable, aeroprediction code was a more 
difllcult task. Nonlinear flow phenomena have normally been studied using higher order 
numerical techniques of the more complex research codes. Semi-empirical computational 
schemes have been described which are applicable to the asymmetric separation prob- 
lem. but many of these approaches are either applied in paneling method codes, or left 
as modeling procedures unincorporated into complete operational missile codes. The 
VORSTAB 11 code may be ideally suited for application on the IBM 3033 or VAX sys- 
tems. but is currently in restricted distribution and not available commercially. Through 
contact with Dr. Lan (University of Kansas Flight Research Laboratory) and Dr. 
.Mehrotra (VlGh'.AN. Inc.), however, the possibility of obtaining a binary version of the 
VORSTAB 11 source code has been discussed. The availability of this code, in any for- 
mat. would provide a greatly enhanced facility to conduct aerodynamic research and 
comparison to experimental data. 

B. SOURCE CODE DESCRIPTION AND OPERATION 

NSWC and MISSILE D.ATCO.M were received on standard 9 track, formatted 
magnetic tapes spools. Each code was accompanied by a users guide, and a theorx' vol- 
ume detailing the various computational methods, data bases, algorithm structure and 
subroutine functions. Test case input and output data were present for both programs. 
The aeroprediction codes were copied from tape to virtual machine disk in order to fa- 
cilitate initial use related to installation and validation; subsequent operations are an- 
ticipated to be performed in a batch mode. Preliniinarx’ work was conducted on the 
Watfor compiler. Compilation in Watfor requires that the source code adhere to much 
more stringent syntax and structure rules than is demanded by VS Fortran compilation. 
While such a cautious approach entails numerous, albeit relatively minor programming 
changes, inadvertant programming errors are greatly reduced. The conversion process 
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additionally allows the operator to become acutely familiar with source code structure 
and subroutine function. This exposure becomes extremely beneficial when trouble- 
shooting is required on larger source codes. 

1. NSWC 

The NSWC program is written in ANSI Fortran, and structured in a top down 
executive manner. Various characteristics and subroutine methodology are presented in 
Table 5, Table 6 on page 22, and Table 7 on page 23; high angle of attack limitations 
are shown in Figure 1 on page 24. 



Table 5. BODY ALONE METHODS OF THE NSWC AEROPREDICTION CODE 





MACH REGION 


COMPONCNT 


SUBSONIC 


TRANSONIC 


LOW 

SUPERSONIC 


HIGH 

SUPERSONIC 


NOSE 

WAVE DRAG 


— 


EULER PLUS 
EMPIRICAL 


SECOND ORDER 
VAN D^'KE PLUS 
MODIFIED 
NEWTONIAN 


SECOND ORDER 
SHOCK 
EXPANSION 
PLUS MODIFIED 
NEWTONIAN 


BOATTAIL 
WAVE DRAG 


— 


WU 

AND 

AOYOMA 


SECOND ORDER 
VAN DYKE 


SECOND ORDER 
SHOCK 
EXPANSION 


SKIN FRICTION DRAG 


VAN DRIEST II 


BASE DRAG 


EMPIRICAL 


INVISCID LIFT AND 
PITCHING MOMENT 


EMPIRICAL 


EULER OR WU 
AND AOYOMA 
PLUS EMPIRICAL 


TSIEN 

FIRST ORDER 
CROSSFLOW 


SECOND ORDER 
SHOCK 
EXPANSION 


VISCOUS LIFT AND 
PITCHING MOMENT 


ALLEN AND PERKINS CROSSFLOW 



Source; [Ref. 36, p. 43) 
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Table 6. WING ALONE AND INTERFERENCE METHODS OF THE NSWC 
AEROPREDICTION CODE 





MACH REGION 


COMPONENT 


SUBSONIC 


TRANSONIC 


LOW' 

SUPERSONIC 


HIGH 

SUPERSONIC 


INVISCID LIFT AND 
PITCHING MOMENT 


LIFTING 

SURFACE 

THEORY 


EMPIRICAL 


LINEAR THEORY 


SHOCK 

EXPANSION 

AND 

STRIP THEORY 


WING-BODY 

INTERFERENCE 


SLENDER BODY 
THEORY AND 
EMPIRICAL 


LINEAR THEORY 
SLENDER BODY 
THEORY AND 
EMPIRICAL 


— 


WING-TAIL 

INTERFLRENCE 


LINE VORTEX 
THEORY 


— 


WAVE DRAG 


— 


EMPIRICAL 


LINEAR THEORY 
AND MODIFIED 
NEWTONIAN 


SHOCK 
EXPANSION- 
STRIP THEORY 
MODIFIED 
NEWTONIAN 


SKIN FRICTION DRAG 


VAN DRIEST 11 


TRAILING FDGE 
SEPARATION WAKE 


EMPIRICAL 


BOD^ BASF 
PR F SSL RE DRAG 
FROM TAIL FINS 


EMPIRICAL 



source; [Ref. 36, p. 4-*] 
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Table 7. DYNAMIC DERIVATIVE COMPUTATION OF THE NSWC 
AEROPREDICTION CODE 





MACH REGION 


COMPO.NENT 


SUBSONIC 


TRANSONIC 


LOW 

SUPERSONIC 


HIGH 

SUPERSONIC 


BODY ALONE PITCH 
DAMPING MOMENT 


EMPIRICAL OR MODIFIED SLENDER BODY THEORY OR 
LINEAR INTERPOLATION OR NEWTONIAN THEORY 


WING AND 

INTERFERENCE ROLl 
DAMPING 


LIFTING 

SURFACE 

THEORY 


EMPIRICAL 


LI.NEAR THEN- 
WING THEORY 


STRIP THEORY 


WING AND 
INTERFERENCE 
MAGNUS MOMENT 


ASSUMED ZERO 


BODY ALONE ROLL 
DAMPING MOMENT 


EMPIRICAL 


WING AND 
INTERFERENCE 
PITCH DAMPING 


SLENDER WING OR NEWTONIAN OR LIFTING SURFACE THEORY 
OR Sl RIP THEORY OR EMPIRICAL 



Source: [Ref. 36, p. 44] 



The source code is approximately 13,500 records in length. Extensive use has been made 
of named COMMON block variables and GO TO statements. This programming style 
allows a marked reduction in the overall source code size, but greatly increases the dif- 
ficulty of troubleshooting and error tracing. The general absence of in-code documenta- 
tion (comments) adds to this dilliculty. The main program body is used to read input 
data, determine subroutine sequencing, and print the majority of output data. Data in- 
put is performed through the creation of a sequential, formatted input file. This method 
is lacking in that examination of the input file provides no direct information as to body 
configuration, flight conditions or run options. Such information is apparent only upon 
comparing the input file with the user manual. Compilation on the IBM system requires 
several CPU seconds, depending on the input configuration and the run options. 

Aside from routine conversion errors, i.e., type specification mismatch, array 
dimensions, initialization, formatting, etc., the NSWC program required ver}' few alter- 
ations to the source code. The majority of changes were connected to index arguments 
in loop structures, CO.M.MON block variable declarations of difierent list element 
length, and precision checks on variable values used in bounded value functions, such 
as square root and hyperbolic functions. Validation of the NSWC program was made 
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(1) Body Geometry'; 

(i) body alone 

(ii) body plus cruciform tail 

(2) Mach Number Range: 0.8 to 3.0 

(3) Angle of Attack: 

(i) isolated component: 0 to 180 degress 

(ii) body with tail: 0 to 45 degrees 

(4) Tail Geometry: 

(i) trapezoidal planform, edges parallel to axis 

(ii) zero sweep trailing edge, parallel to base 

(iii) leading edge sweep from 0 to 70 degrees 

(iv) taper ratio 0 to 1.0 

(v) aspect ratio (two fin) 0.5 to 2.0 

(5) Nose Length: 1.5 to 3.5 calibers (pointed ogive) 

(6) Afterbody Length: 6 to 18 calibers 

(7) Total Span-to-Diameter Ratio; 1.0 to 3.33 

Figure 1. High Angle of Attack Limits for the NSWC Aeroprediction Code 

against 1 1 test case input output data sets. These inputs represent a wide variety of 
standard configuration models used in aerodynamic research. The test case input is de- 
signed to cover the various input options allowable, and thus perform an operability 
check on nearly every subroutine. Of the 1 1 trial runs, the 9 involving ogive or conical 
nose sections were completed successfully, with good accuracy for all geometry param- 
eters. structural loading values and aerodynamic quantities listed in the test case output. 
The blunt or truncated nose configuration output could not be replicated. For the given 
geometry and reference conditions, the pressure coefficients (with respect to axial body 
position) were seen to fluctuate in a seemingly random fashion. Discussions with Dr. 
Devan indicate that the problem may be related to the (double precision) accuracy of the 
solution as calculations are made to determine the second order Van Dyke jump condi- 
tions across slope and curvature discontinuities [Ref 53]. It is felt that this problem can 
be corrected. 
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2. MISSILE DATCOM 

The MISSILE DATCOM source code is written in ANSI Fortran, and has the 
appropriate syntax included to conform to the coding standards of both Fortran IV and 
Fortran V (Fortran 66 and Fortran 77) compilers. At roughly 42,000 records, the source 
code size represents the upper limit for manageable operation on the interactive 
timeshare net. Storage requires a single designated disk; as the source code is too large 
to edit as a single body, the program was subsequently divided into two separate files 
of equal length. Compilation involves copying both files to a temporary' access disk, 
appending the files into the original source code, and executing the program. If changes 
must be made to the source code, the process must be repeated. The Watfor compiler, 
with the Release 31 option, is able to execute the program, although 4 .MB and 20 cyl- 
inders of temporar}' access disk are required. During high volume loading of the inter- 
active timesharing net, a routine compilation might be expected to exceed 25 minutes 
of real time. Initial installation and validation procedures were performed during low 
system usage periods. 

MISSILE DATCOM is structured in a top down executive manner. Body alone 
methods are shown in Table 8 on page 26. Fin alone and configuration synthesis tech- 
niques are presented in Table 9 on page 27. 
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Table 8. BODY ALONE METHODS OF THE MISSILE DATCOM CODE 





MACH REGION 


COMPONENT 


SUBSONIC 


TRANSONIC 


SUPERSONIC 


NOSE 

WAVE DRAG 


— 


EULER PLUS 
EMPIRICAL 


VAN DYKE HYBRID 
PLUS MODIFIED 
SECOND ORDER 
SHOCK EXPANSION 


BOATTAIL 
WAVE DRAG 


— 


PAYNE 

CORRELATION 


VAN DYKE HYBRID 
PLUS MODIFIED 
SECOND ORDER 
SHOCK EXPANSION 


SKIN 

FRICTION DRAG 


BLASIUS PLUS TRANSITIO.N’ PLUS VA.\ DRIEST II 


BASE DRAG 


MODIFIED SECOND ORDER SHOCK EXPANSION 
PLUS EMPIRICAL 


INVISCID LIFT AND 
PITCHING MOMENT 


EMPIRICAL 


EMPIRICAL 


VAN DYKE HYBRID 
SECOND ORDER 
SHOCK EXPANSION 


VISCOUS LIFT AND 
PITCHING MOMENT 


BAKER PLUS JORGENSEN CROSSFLOW 



Source; (Ref. 26, p. IS) 
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Table 9. WING ALONE AND INTERFERENCE METHODS OF THE MISSILE 
DATCOM CODE 





MACH REGION 


COMPONENT 


SUBSONIC 


TRANSONIC 


SUPERSONIC 


INVISCID LIFT AND 
PITCHING MOMENT 


LOWRY-POLHAMUS 

CORRELATION 


R.A.S. DATA SHEET 
INTERPOLATION 


SUPERSONIC 
WING THEORY 


CARRY-OVER 

INTERFERENCE 


NACA 1307 


VIRA AND FAN 
CLOSED FORM 
ANALYTIC SOLUTIONS 


WAVE DRAG 


— 


LINEAR FAIRING 
(FROM M = 1.05) 


LINEAR THEORY 
SHOCK EXPANSION 
STRIP THEORY 
MODIFIED NE\\TONTA^ 


SKIN FRICTION DRAG 


BLASIL'S PLUS TRANSITION PLUS VAN DRIEST II 


'TRAILING EDGE 
SEPARATION 


EMPIRICAL 


BODY BASE 
PRESSl RE DRAG 
DLE TO TAIL FINS 


E.MPIRICAL 



Source: [Ref. 26, p. 70] 



Data input is accomplished using a namelist construct. Although the namelist format is 
dilTerent than that required of Fortran compilers, a namelist emulator is incorporated 
into the code. This namelist input method is simple to use and allows a much better case 
documentation than sequential data input files. Control cards are available to specify 
which options are desired for the compilation. The MISSILE DATCO.M test case pro- 
vides input output data for an arbitrary body- wing-tail missile. While this test case al- 
lows verification of certain principal subroutines, the complexity of this code is such that 
a more complete test series would prove beneficial. The majority of errors encountered 
during the validation procedure were related to variable type declaration, array dimen- 
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sions, and COMMON block equivalence. The program has an internal error checking 
feature which is of great assistance in troubleshooting. The extensive use of documen- 
tation (comment) statements is excellent, and nearly essential in the maintenance of any 
large code. Only one structural change was necessary to facilitate compilation of the 
program. This alteration removed a logical condition from a subroutine call statement. 
The condition was deleted in order to allow variable loading and array initialization 
which were otherwise undefined. 

The test case output data was matched with exceptionally good accuracy. This 
is due, in part, to the source code design which declares the majority of variables and 
subroutines as implicit double precision for non-CDC machines. One source of con- 
fusion is present in the user guide. As the "DA.MP" option was not a control card input, 
dynamic derivative values are not listed within the test case output. Dynamic derivative 
output is shown, however, as a separate page to illustrate the format type. The reference 
conditions listed for this format example strongly suggest that the values reflect the 
output computation based on the test case input data, although these sums could not 
be replicated. Further ambiguity in the prediction of dynamic derivatives is discussed in 
the data comparison section of this paper; prediction methods are displayed in 
Table 10 on page 29. 
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Table 10. DYNAMIC DERIVATIVE METHODS OF THE MISSILE DATCOM 
CODE 





METHOD 


DERIVATIVE 


BODY 


FIN 


CNq 


USAF STABILITY AND 
CONTROL DATCOM 


USAF STABILITY AND 
CONTROL DATCOM 


Cmq 


USAF STABILITY AND 
CO.NTROL DATCOM 


USAF STABILITY AND 
CONTROL DATCO.M 


CNAD 


USAF STABILITY AND 
CONTROL DATCOM 


USAF STABILITY AND 
CONTROL DATCOM 


Cm AD 




USAF STABILITY AND 
CONTROL DATCOM 


Cmq + CmAD 


ERICSSON (LMSC) 


Cmq ^ CmAD 


Op 


ARDC SPINNER CODE (SPIN-73) 




Cnp/sin(alpha) 


ARDC SPINNER CODE (SPIN-73) 




CIp 


ARDC SPINNER CODE (SPlN-73) 




Source: (Ref. 26, p. 106 
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IV. COMPARISON OF PREDICTION AND RESULTS 



Following installation and successful test case validation procedures, each 
aeroprediction code was prepared for comparison runs. In order to permit a meaningful 
interpretation of source code performance, an acceptable experimental data base was 
required as a benchmark for each trial run; comparison between the MISSILE 
DATCOM and NSWC programs was dependent on a common description of body ge- 
ometr>' under equivalent reference conditions. 

Trial run inputs were selected from the test case validation section of the NSWC 
program users guide [Ref 54]. This input data covers a sufficiently wide variety of missile 
and projectile configurations, and is accompanied by experimental wind tunnel and 
range data for the appropriate flight conditions. Wind tunnel data and geometry inputs 
of the .Aeronautics and Astronautics Department Standard-type missile model were also 
available. The following sections provide a comparison of NSWC and MISSILE 
DATCOM source code predictions of various aerodynamic coefficients. 

A. ARMY-NAVY SPINNER 

The Army-Na\w Spinner is a spin stabilized projectile geometiy’ which is illustrated 
in Figure 2. The comparison trials were made at zero angle of attack with no roll angle 
and no control dellections; .Mach condition was incremented from subsonic to high 
supersonic. Geometry inputs to both programs were identical, excepting rotating band 
height and spin stabilization, which can not be included into the MISSILE D.ATCOM 
input file. The comparison experimental data appearing in Figure 3 through Figure 6 are 
from Arnold Engineering and Development Center (,AEDC) (Ref 55]. 

1. Normal Force Coefficient Derivative 

The prediction values of the normal force coefficient slopes are presented in 
Figure 3 on page 32. The NSWC curve is seen to follow the .AEDC data fairly well, al- 
though underprediction and overresponse are present in the transonic region. Prediction 
accuracy improves with higher .Mach numbers. 

The .MISSILE DATCO.VI prediction curve reflects the general form of the 
AEDC data and the NSWC curve, but shows a continuous underprediction in all .Mach 
regions. Prediction within the transonic region is a better fit, in magnitude, than the 
NSWC code, but appears to lag slightly behind as a function of .Mach number. The 
rapid increase in slope at the transonic-supersonic juncture is only partially present. 
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The discrepancy in prediction between these codes is most probably a result of 
the different methods used to calculate normal force and pitch damping coefficients. 
Within the transonic region, the NSWC code employs an improved method of Nielsen 
Engineering and Research (NEAR) [Ref 56]. This method, which is restricted to con- 
stant radius afterbodies of five caliber length of less, is based on (one degree angle of 
attack) normal force and pitch moment solution data of the Euler equations. Using a 
least squares functional fit to a truncated Taylor series expansion, pressure distribution 
is interpolated from the various configuration grid data as determined by NE.-\R. The 
original interpolation scheme has been improved by augmenting the functional form fit 
data with computations from an improved Van Dyke Hybrid potential model for the 
transonic-supersonic crossover point. This additional data were generated from a much 
larger data base, and significantly improves the accuracy, since not all of the NE.AR 
Euler solutions were completely convergent. A less accurate method of Moore is retained 
for use with afterbodies exceeding five calibers [Ref 57]. [Ref 36, pp. 32-43, 64). The 
action of this computational method is clearly seen at the Mach 1.2 portion of the 
NSWC curve in Figure 3. 
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Figure 3. Normal Force Cune Slope Comparison for (he Army-Na\y Spinner 
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The MISSILE DATCOM program uses empirical methods to estimate the 
normal force and pitch moment coefficients in the subsonic and transonic regions; 
modified second order shock expansion and Van Dyke Hybrid techniques are incorpo- 
rated for supersonic Mach. The NEAR method of Ref. 56 was considered during pre- 
liminary feasibility recommendations, but was rejected in view of the wider range of 
fineness ratios available from the empirical Messersmitt-Bolkow-Blohm (.MBB) data 
sheets. The MISSILE DATCOM algorithms are reportedly far superior to those found 
in the NSWC code; however, this claim is not supported by the trial run comparison of 
the Army-Na\y Spinner. (Ref 26, pp. 19-20, 30-32]. 

2. Center of Pressure 

NSWC and MISSILE DATCO.VI prediction curves for the center of pressure 
are shown in Figure 4 on page 34. The prediction data are plotted as a ratio to body 
length, and referenced 2.51 calibers from the nose tip. The prediction curves are similar 
in form to those in Figure 3, which is expected in light of the previous discussion. Pre- 
diction performance of both codes is somewhat poor for subsonic and low transonic 
Mach number. Significant underprediction error is again evident in the MISSILE 
DATCO.M predictions within the supersonic region. 

3. Drag Coefficient 

The zero angle of attack drag predictions are presented in Figure 5 on page 
35. As can be seen, both NSWC and MISSILE DATCOM show close correlation to 
the AEDC experimental data. This level of agreement is not unexpected, however, in 
that both programs use similar computational methods for nose and body components. 
Transonic prediction is performed by an improved wave drag prediction model of NEAR 
[Ref 5Sj. This technique estimates drag efl'ects through interpolation of pressure dis- 
tribution data generated by time asymptotic solutions of the Euler equations. The data 
set of pressure distributions was developed from a family of tangent ogives for transonic 
Mach numbers. Short nose and small bluntness computations have been improved by 
recomputing wave drag with the ICANBOD program [Ref 59]. Supersonic prediction is 
accomplished through a potential flow or a second order shock expansion method. The 
second order Van Dyke and modified second order shock expansion techniques of the 
NSWC code provide a better fit for supersonic .Vlach predictions, with drag coefficient 
estimates nearly coincident with the experimental data. [Ref 36, pp. 27-29, Ref 26, pp. 
23-26]. 
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Figure 4. Center of Pressure Comparison for the Army-Na >7 Spinner 
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(ALPHA=0) 

Figure 5. Drag Coefficient Comparison for the Army-Nan’ Spinner 
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4. Pitch Damping Coefncient 

The NSWC, MISSILE DATCOM and AEDC pitch damping coefficient data 
appear in Figure 6 on page 37. The MISSILE DATCOM prediction values have been 
increased by a factor of two to account for differences in the non-dimensional definition 
of the pitch damping coefficient. Predicted values of the two codes are nearly identical 
in the subsonic region, although the fit is marginal in comparison to the AEDC data. 
Low supersonic prediction of the NSWC code is superior to that of MISSILE 
DATCO.M; however, neither program offers a particularly accurate estimate within the 
supersonic region. The divergence of the two prediction code curves may result from the 
different computational methodology. 

Dynamic derivative prediction of the NSWC code is strictly valid for low angles 
of attack. The prediction technique is a modified application of the L.MSC method of 
Ericsson (Ref 60]. The LMSC method has shown good accuracy for subsonic and 
transonic regions, where slender body theory remains applicable. Accuracy is diminished 
at supersonic Mach numbers, however, since boattail and afterbody effects are not 
considered. In order to reduce the prediction error within the supersonic region, the po- 
tential flow construct is assumed invalid, and solution is supplemented by strip theory. 
While this modified LMSC approach is adequate for most conditions, the prediction 
method is further augmented by an empirical process derived from wind tunnel and 
ballistic range test data on spin stabilized projectiles. This program is known as SPIN- 
NER or SPlN-73 (Ref 61 j. For body alone configurations, the NSWC code is structured 
to compare prediction values of the LMSC and SPINNER methods; if the LMSC pre- 
dicted pitch damping coefficient is less negative and not within 75 percent of the SPIN- 
NER value, the SPINNER value is chosen. The effect of this empirical method results 
in the significant difference between the NSWC and .MISSILE D.\TCOM prediction 
values at supersonic .Mach numbers. [Ref 36, pp. 21-25. 5S|. 

The .MISSILE D.ATCO.M code uses a follow-on version of the LMSC method, 
but does not incorporate any compensating modifications for the prediction of pitch 
damping coefficients [Ref 62]. The SPINNER routine is present, but contributes only 
to the solution of yawing moment, rolling moment and side force due to roll rate, as 
shown in Table 10 on page 29. 

B. BASIC TINNER 

The Basic Finner model is a cone-cylinder body with tail configuration as presented 
in Figure 7. Comparison computer runs were conducted at zero degree angle of attack 
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Figure 6. Pitch Damping Coefficient Comparison for the Army-Na\y Spinner 
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with no control deflection and no roll angle. Input conditions included subsonic, 
transonic and supersonic Mach numbers. The experimental data in Figure 8 through 
Figure 11 are from a Bureau of Weapons technical report (Ref. 63]. 

1. Normal Force Coefficient Derivative 

The XSWC and MISSILE DATCOM prediction values are plotted against the 
experimental data in Figure 8 on page 40. The prediction curves of both programs show 
a fairly accurate fit with the experimental data. The MISSILE DATCOM predictions 
are noticeably superior to the XSWC values, and are nearly identical to the experimental 
results within the supersonic region. The XSWC curve reflects a roughly constant 
underprediction error for all Mach numbers. While both codes show a greatly improved 
prediction capability relative to the Army-X'a \7 Spinner projectile of Figure 3, it is ap- 
parent that the MBB data used by MISSILE DATCOM are particularly better suited 
for application to more traditional missile geometries. The XSWC code shows some 
underprediction error for both the Spinner and Basic Finner configurations. 

2. Center of Pressure 

The center of pressure comparison is shown in Figure 9 on page 41. Coefficient 
values are expressed as ratios to body length, referenced 6. 10 calibers from the nose. The 
prediction accuracy of the XSWC and MISSILE D.ATCOM codes is excellent, especially 
within the low supersonic Mach region. As with the normal force derivative predictions, 
the estimation of center of pressure coefficients shows dramatic improvement from the 
Army-Xavy Spinner case. The increased accuracy is most pronounced for the MISSILE 
D.ATCOM output. Both programs reflect an underprediction trend within the high 
supersonic region, although to a much smaller degree than for the Army-Xavw Spinner 
comparison of Figure 4. The XSWC code provides a better functional fit for transonic 
Mach numbers. 

3. Drag Coefficient 

Prediction curves and experimental data of zero degree angle of attack drag co- 
efficients are presented in Figure 10 on page 43. Both programs are seen to underpre- 
dict drag coefficient magnitude within the transonic Mach region, although the XSWC 
curve provides a better fit than does the MISSILE DATCO.M curve. The prediction 
accuracy of each code improves for supersonic .Mach numbers, with the XSWC predic- 
tion curve maintaining a better correlation to the experimental data. The quality of drag 
coefficient prediction is significantly reduced as compared to the Army-Xa\y Spinner 
projectile of Figure 5. As previously discussed, the XSWC and .MISSILE D.ATCOM 
codes employ nearly identical computational methods for nose and cylinder drag calcu- 
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lation. While a reduction in prediction accuracy might be expected as the complexity 
of the geometry increases, it is reasonable to assume that the reduced similarity between 
the NSWC and MISSILE DATCOM drag predictions results from the calculation 
methods for wing alone drag effects. Both prediction programs utilize identical empirical 
methods for fin alone computations of trailing edge separation drag and skin friction 
drag. A difference exists, however, in the calculation of wave drag effects. The MISSILE 
D.ATCOM code employs the methodology of Moore in both the transonic and super- 
sonic .Mach regions [Ref 64]. Transonic wave drag is estimated through a potential flow 
solution at Mach number 1.05, which is then linearly reduced as a function of Mach 
number. The supersonic technique is a computational grid solution using potential flow 
theory-. [Ref 26, pp. 70-75). 

The NSWC program uses these same methods of .Moore, but only within the 
supersonic Mach region; transonic wave drag effects are estimated using an empirical 
method. From the comparison of drag prediction in Figure 10, it is concluded that the 
empirical method of the NSWC code is superior to the approximate potential flow (lin- 
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Figure 8. Normal Force Cune Slope Comparison for the Basic Finner 
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Figure 9. Center of Pressure Comparison for the Basic Finner 
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ear fairing) method of the MISSILE DATCOM code for drag coefTicient estimation at 
transonic Mach numbers. 

4. Pitch Damping Coefficient 

Pitch damping coefficient predictions are presented in Figure 1 1 on page 44. 
While the NSWC curve reflects a fairly accurate image of the experimental data, it is 
obvious that the MISSILE DATCOM values are seriously underpredictive. Upon com- 
pilation of the .MISSILE DATCO.M comparison case for Basic Finner, an examination 
of the source code was conducted in an effort to uncover possible structure errors. The 
inclusion of an installation related error was considered a possibility, especially since no 
dynamic derivative calculations are included in the MISSILE DATCOM validation test 
case. The inspection of the source code failed to reveal any obvious errors or inadvertent 
alterations. 

In reviewing the .MISSILE DATCO.M Final Report (Ref 26). a figure contain- 
ing pitch damping prediction data for the Basic Finner was located. The model geometry 
is also shown, and is identical to the input configuration for the comparison runs! 
MISSILE DATCOM prediction values and experimental data are plotted in this figure 
for two different moment reference positions; however, the prediction values could not 
be replicated with subsequent compilations of the .MISSILE DATCOM source code. 
Furthermore, the unreferenced experimental data are significantly different from the data 
of Ref 63 which appears in the NSWC User Manual. Finally, the MISSILE DATCOM 
Basic Finner comparison values presented here in Figure 1 1 appear similar to neither the 
experimental data nor program output curves shown in the figure on page 63 of the 
.MISSILE DATCO.M Final Report. 

A further attempt to resolve these discrepancies was made using a dynamic de- 
rivative output sheet on page 123 of the .MISSILE DATCO.M User .Manual (Ref 35). 
Although these data are not part of the validation test case, an examination of the flight 
conditions and reference quantities for these data strongly suggests that the dynamic 
derivative values reflect prediction output for the test case input geometiy. The test case 
was rerun accordingly, with the dynamic derivative control card option; however, the 
User .Manual values could not be generated. It is somewhat interesting to note that the 
magnitude of error for this test case run is roughly equivalent to that e.xperienced for the 
Basic Finner comparison. While this problem requires further attention, the validity of 
the experimental Basic Finner data (Ref 63) used by the NSWC program has been 
confirmed through alternate sources [Refs. 65,66]. No interpretation or comparison of 
the .MISSILE DATCO.M pitch damping coefficient prediction data is offered. The 
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MISSILE DATCOM User Manual pitch damping output and the Final Report figure 
of Basic Finner pitch damping coefficient prediction are enclosed in Appendix A and 
Appendix B. 

C. AIR SLEW DEMONSTRATOR 

The Air Slew Demonstrator is a tangent ogive body with tail as illustrated in 
Figure 12 on page 46. The comparison runs were made at low supersonic speed 
(M= 1.3) for angles of attack between 0 and 50 degrees. Control deflection and roll 
angle were not considered for normal force and center of pressure predictions; roll mo- 
ment was predicted for a roll angle of 22.5 degrees. The experimental data shown in 
Figure 13 through Figure 15 are from AEDC [Ref 67]. 

1. Normal Force Coefficient 

Normal force coefficient predictions are presented in Figure 13 on page 47. 
Both the NSWC and MISSILE DATCOM prediction curves show a good degree of 
accuracy relative to the experimental data. The NSWC prediction values appear to re- 
flect a roughly constant underprediction; the MISSILE DATCO.M curve shows an 
overprediction trend which is largest between 25 and 45 degrees. The NSWC code uses 
a Martin Marietta empirical routine for the estimation of high angle of attack effects 
[Ref 6Sj. The high angle of attack data set is maintained as a separate input file or tape. 
.VIISSILE DATCOM incorporates the Allen and Perkins plus Jorgensen crossflow 
method to approximate high angle of attack effects. This technique provides generally 
accurate predictions, with a maximum deviation from experimental data in the 15 to 60 
degree angle of attack range [Ref 26, pp. 30-31.]. The comparison of normal force co- 
efficients for the Air Slew Demonstrator indicates that the NSWC and MISSILE 
DATCOM programs are quite comparable for transonic, high angle of attack condi- 
tions. The NSWC code appears to provide a marginally better estimate of the nonlinear 
trend at angles of attack greater than 25 degrees, which is an incidence angle normally 
associated with the onset of asymmetric vortex shedding. 

2. Center of Pressure 

Center of pressure coefficient predictions appear with the AEDC experimental 
data in Figure 14 on page 48. These data are expressed as ratios to the body length, 
referenced 4.33 calibers from the nose. Unlike the zero degree angle of attack compar- 
isons, the high angle of attack predictions of center of pressure are, at best, considered 
fair to marginal. The NSWC curve fails to represent the experimental data for angles of 
attack less than 30 degrees; higher angle of attack estimates are overpredicted but of a 
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Figure 12. Air Slen Demonstrator (Ref. 54, p. 60] 

functional form similar to that of the experimental data. While the MISSILE DATCOM 
curve is overpredictive for all angles of attack, the accuracy is improved below 15 degrees 
angle of attack. The parabolic form of the NSWC curve between 5 and 25 degrees is 
significantly in error relative to the experimental data and the .MISSILE DATCOM 
predictions. 

3. Roll .Moment Coefficient 

Roll moment calculations were made with a 22.5 degree roll angle input. The 
NSWC and .MISSILE D.ATCO.M predictions are presented with experimental data in 
Figure 15 on page 50. The NSWC program uses an empirical data base for roll moment 
prediction, while .MISSILE DATCOM employs the SPIN-73 program. Roll moment 
estimations are of the right order of magnitude for both programs, but neither code is 
particularly accurate. The .MISSILE DATCO.M curve shows a significant overprediction 
trend as a function of increasing angle of attack. The NSWC curve shows a superior 
prediction capability at higher angle of attack, which is indicative that the empirical data 
base of the NSWC code provides better accuracy than the SPIN-73 routine at high an- 



46 



NORMAL FORCE COEFFICIENT 

0.0 2.0 4.0 6.0 0.0 10.0 12.0 14.0 16.0 



□ 


LEGEND 

MISSILE DATCOM 


o 


NSWC 


A 


EXPERIMENTAL DATA 




[Ref. 67 j 




Figure 13. Normal Force Coefficient Comparison for the Air Slew Demonstrator 
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Figure 14. Center of Pressure Comparison for the Air Sle>\ Demonstrator 
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gles of attack. A rough equivalence exists between the two programs for prediction at 
angle of attack less than 40 degrees. 

D. TMX-2774 (T-9 TAIL) 

The TMX-2774 geometr}’ is presented in Figure 16 on page 51. Comparison testing 
was conducted at zero degree angle of attack, with no control deflection and no roll 
angle. Mach conditions were incremented from low supersonic to high supersonic. The 
experimental data in Figure 17 and Figure 18 are from Nichols [Ref. 69]. 

1. Normal Force Coefficient Derivative 

Normal force curve slope predictions for the TMX-2774 appear in Figure 17. 
The NSWC curve shows slight underprediction for the low supersonic region and a slight 
overprediction for higher Mach numbers. The MISSILE DATCOM curve maintains an 
overprediction for all Mach numbers, although the low supersonic Mach estimates are 
a much closer fit to the experimental data. The overall accuracy of the NSWC code is 
good, while the MISSILE DATCOM program can be considered less correct at high 
Mach numbers. 

2. Center of Pressure 

The center of pressure predictions are plotted with the experimental data in 
Figure 18 on page 53. The data are taken as ratios to the body length, and referenced 
6.81 calibers from the nose. The NSWC prediction curve is of excellent accuracy with 
respect to the experimental data. The MISSILE DATCOM estimates are of the right 
functional form, but are underpredicted for all Mach conditions. 

E. TMX-I751 

The TMX-1751 is a wing-body-tail configuration, as is shown in Figure 19 on page 
54. The comparison runs were conducted within the supersonic Mach region at zero 
degree angle of attack. Control dellections and roll angles were not considered. Exper- 
imental data for Figure 20 through Figure 22 are from Ref 69. 

1. Normal Force Coefficient Derivative 

The normal force curve slopes are presented in Figure 20 on page 55. The 
NSWC and MISSILE DATCOM values show good correlation to the experimental 
data, particularly between Mach 2.5 and .Mach 4.5. The NSWC program provides a high 
level of accuracy within this Mach range; however, a small degree of overprediction is 
apparent. The MISSILE DATCOM prediction curve is also accurate within this Mach 
range, although the amount of overprediction is noticeably greater than for the NSWC 
program. 



49 



n 


LEGEND 

MISSILE DATCOM 


o 


NSWC 


A 


EXPERIMENTAL DATA 




(Ref. 67] 



o 




M=1.30 PH1=22.5 



Figure 15. Roll Moment Comparison for the Air Slew Demonstrator 
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2. Center of Pressure 

The predictions of center of pressure coefficients are presented in Figure 21 on 
page 56. The output values are expressed as ratio to the body length, referenced S.41 
calibers from the nose tip. The MISSILE DATCO.M and NSWC plots are nearly iden- 
tical in form, with both programs underpredicting the experimental data. While the 
.MISSILE DATCO.M and NSWC programs fail to replicate the functional form of the 
experimental data, the prediction values are of the correct order of magnitude, such that 
the relative errors are reasonably small. In comparison to the trend of the experimental 
data, the linear form of the NSWC and MISSILE DATCOM curves suggest that the low 
supersonic prediction accuracy would be increasingly poor. 

3. Drag Coefficient 

The zero degree angle of attack drag coefficient data are shown in Figure 22 
on page 57. Both the NSWC and MISSILE DATCOM programs overpredict drag co- 
efficient values in comparison to the experimental data, although the functional form 
of the prediction curves is basically correct. The degree of error is greater than in the 
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Figure 17. Normal Force Curve Slope Comparison for the TMX-2774 (T-9 Tail) 
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Figure 18. Center of Pressure Comparison for the TMX-2774 (T-9 Tail) 
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Figure 19. TMX-1751 [Ref. 54. p. 67) 

body-tail Basic Finncr comparison, with the trend to overpredict tlie drag coelFicents 
rather than the underprediction seen for the Basic Finner. 

F. STANDARD MISSILE MODEL 

The Standaid Missile model configuration is presented in Figure 23 on page 58. 
1 his scale model is sintilar in design to the Navy Standard Missile, and is used in support 
of aerod\namic research at the Naval Postgraduate School. Previous work by Rabang 
included high angle of attack wind tunnel testing of this model in wing-body-tail and 
body-tail configurations; the empirical data of this testing are used as the reference data 
in the following comparisons [Ref 70]. 

Comparison runs were made at a low subsonic .Mach number of 0.10 and a 
Reynolds number of 1 . 1 E5. Angle of attack was incremented from 0 to 90 degrees. These 
conditions prevented use of the NSWC code, which is restricted to transonic and super- 
sonic Mach for high angle of attack analysis, as shown in Figure 1 on page 24. Body- 
wing-tail and body-tail configurations were input with no roll angle; the body-wing-tail 
configuration was additionally run with a 45 degree roll angle input. 
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Figure 20. Normal Force Curve Slope Comparison for the TMX-1751 
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Figure 21. Center of Pressure Comparison for the TMX-1751 
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Figure 22. Drag Coefficient Comparison for the TMX-1751 
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NAVAL POSTGRADUATE SCHOOL 




Geometn inputs for tlie body-tail configuration case were rather direct; however, 
the description of the wing geornctrv' was more difficult due to the non-traditional shape. 
An appro.ximate fin description was input which ignored the break in the wing leading 
edge, 1 his was done by defining the tip chord dimension at the tip span station such that 
a continuous leading edge is taken from the leading station of the root chord at a sweep 
angle of 30 degrees. 

1. Normal Force Coefficient for the Body-lail Configuration 

1 he MISSILE DATCOM prediction is quite good, as can be seen in Figure 24 
on page 59. The normal force coefficient is very accurately described up to 50 degrees, 
at which point an underprediction trend commences. As the slope of the experimental 
data curve begins to rapidly diminish, near 65 degrees, the MISSILE DATCOM curve 
maintains a strong positive slope. The high angle of attack portion of the MISSILE 
DATCOM curve appears functionally correct, but somewhat unresponsive to the non- 
linear effects reflected in the wind tunnel data for angles of attack greater than 50 de- 
grees. This angle of incidence corresponds to the observed onset of unsteady, asymmetric 
separation. 
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Figure 24. Normal Force Comparison for the Standard Missile Model 
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2. Normal Force Coefficient for the Body- Wing-Tail Configuration 

As can be seen from the data curves in Figure 25 on page 61, the MISSILE 
DATCOM estimates of normal force coefficient are underpredicted for the more com- 
plex body-wing-tail configuration. The prediction error is maximum between 15 and 65 
degrees, which is a knoun range of deviation for the Allen and Perkins plus Jorgensen 
crossflow method, as was previously mentioned for the Air Slew Demonstrator. This 
angle of attack range is roughly equivalent to the incidence envelope in which nonlinear 
aerodynamic effects are most pronounced. It is unclear what effect the approximate wing 
description contributed to the relative prediction error, although it can be seen that the 
inclusion of wing effects dramatically shifts the onset of nonlinear lift to lower angles 
of attack. 

3. Normal Force Coefficient for the Body- Wing-Tail Configuration 

The normal force curves shown in Figure 26 on page 62 indicate a nearly con- 
stant underprediction for the Standard .Missile model at a 45 degree roll angle. .MISSILE 
D.ATCOM accuracy is fairly good up to 70 degrees, above which the prediction curve 
continues at a positive slope without the inflection and down turn present in the wind 
tunnel data. 

G. SUMMARY OF RESULTS 

'fhe comparison runs between the NSWC and .MISSILE DATCO.M aeroprediction 
codes indicate that both programs are suitable for preliminary design applications in 
which an accurate prediction of the trend in aerodynamic coefficients (at the correct or- 
der of magnitude) is of greater importance than the precision of specific prediction val- 
ues. The following sections generalize the prediction performance of each code, based 
on an examination of the output from the proceeding comparison trials. 

1. MISSILE DATCOM 

The .MISSILE D.ATCO.M prediction of normal force coefficients is considered 
reasonably accurate for both low and high angle of attack conditions. The best predic- 
tion was output for the Basic Finner model in Figure 8 on page 40, although the normal 
force estimates for the Air Slew Demonstrator. T.M.X-1751 and Standard .Missile model 
are nearly as accurate. The normal force coefficients for high angle of attack conditions 
are seen to be slightly overpredicted for the Air Slew Demonstrator and Standard .Missile 
model (body-tail) configuration. A significant underprediction is present for the body- 
wing-tail Standard .Missile model at angles of attack normally associated with nonlinear 
aerodynamics. 
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Figure 25. Normal Force Comparison for the Standard Missile Model 
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Figure 26. Normal Force Comparison for the Standard Missile Model 
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Center of pressure predictions are accurate to within 10 percent of the exper- 
imental data for the TMX-2774 and Basic Finner, while the Army-Navy Spinner and 
Air Slew Demonstrator estimates show prediction errors of up to 100 percent. An 
underprediction error is characteristic, particularly for high supersonic Mach numbers. 
The estimates of axial position are, however, generally of the correct order of magnitude, 
and typically reflect the functional form of the experimental data. The most accurate 
prediction of center of pressure was for the Basic Finner, while the high angle of attack 
predictions for the Air Slew Demonstrator show significant error. 

The most accurate drag coefficient predictions are within a few percent of the 
experimental data, while the less correct values reflect an 85 percent error. Accuracy is 
diminished within the transonic and supersonic regions. Above Mach 2.5, the drag co- 
efficient estimates are seen to converge on the experimental data. The quality of the drag 
coefficient predictions is seen to decrease as the configuration geometry becomes more 
complex. The Army-Na \7 Spinner comparison of Figure 5 on page 35 shows a good 
level of accuracy with the experimental data. Drag coefficient estimates for the cone- 
cylinder-tail configuration of the Basic Finner are much less accurate, but still provide 
a loose correlation to the empirical values, as shown in Figure 10 on page 43. The 
body-wing-tail geometrx’ of the TMX-1751 results in still further reduction in the accu- 
racy of the drag coefficient predictions, with the low supersonic portion of the output 
divergent from the experimental data. 

An evaluation as to the performance of dynamic derivative prediction is not 
possible. The pitch damping coefficient data in Figure 6 on page 37 and Figure 11 on 
page 44 show excessive underprediction of the experimental data, although a trend sim- 
ilar to the NSWe curve and the reference data is displayed. As there is some degree of 
ambiguity surrounding the MISSILE DATCOM dynamic derivative output, and no test 
case or experimental data which can be used for source code verification or comparison, 
it must suffice to report that the program currently provides no reliable prediction of 
pitch damping coefficients. It is believed, however, that the computational methods are 
theoretically sound, and that the error is connected to the reference definition of the 
output quantities. 

2 . NSWC 

Normal force coelTicient prediction is considered ver>’ accurate. Prediction val- 
ues are within 10 percent of the experimental data, and consistently replicate the func- 
tional form of the empirical results. Some underprediction is present in the normal force 
coefficient estimates of each configuration trial; however, the degree of error is charac- 
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teristically small. No significant degradation in the quality of prediction is seen to result 
from model complexity or angle of attack; a modest increase in prediction error occurs 
for transonic Mach numbers in the Basic Finner output of Figure 8 on page 40. as well 
as for the Army-NaNw Spinner of Figure 3 on page 32. The high angle of attack pred- 
ictions for the Air Slew Demonstrator are exceptionally accurate. 

Center of pressure is predicted with less accuracy than the normal force coeffi- 
cient. While the Basic Finner and T.MX-2774 comparisons fit reasonably w'ell with the 
experimental data, the Air Slew Demonstrator and TMX-1751 predictions are accurate 
only in magnitude. Center of pressure coefficients are generally underpredicted. The ac- 
curacy of prediction at angle of attack is marginal, as indicated by the Air Slew’ Dem- 
onstrator; Figure 14 on page 48 projects a confusing movement for the center of 
pressure at angles of attack between 5 and 30 degrees which can not be adequately ex- 
plained. 

Drag coefficient predictions do not consistently replicate the experimental data. 
The prediction error is seen to vary from just a few percent for the Army-Na\'y Spinner 
projectile, to roughly 80 percent for the TMX-1751. The accuracy is seen to diminish as 
the configuration geometry of the model becomes more complex. While accuracy is ex- 
cellent for the Army-Xa\w Spinner projectile, an underprediction error is evident for the 
body-tail Basic Finner; overprediction of drag coefficients results for the body-wing- tail 
TMX-1751. In all cases, however, the functional form of the drag coefficient prediction 
curves is in close agreement with the experimental data. 

Pitch damping coefficient estimates are satisfactory, although only two config- 
uration cases were examined. The coefficient values are typically overpredicted for all 
Mach regions, although the trend of the prediction curve is quite similar to the exper- 
imental data. The pitch damping coefficient predictions for the Basic Finner are within 
15 percent of the experimental data; a 40 percent error is present for the high Mach 
predictions of the Army-Natw Spinner. Subsonic and transonic predictions are more 
accurate than those for the supersonic Mach region. The predicted values for the Basic 
Finner are in very close agreement with the experimental data; the overprediction trend 
at high Mach numbers for the .Army-NaNw Spinner may well be the result of the com- 
parative value selection process of the LMSC and SPINNER routines. 
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V. CONCLUSIONS 



The usefulness or worth of aeroprediction codes can be measured in various ways. 
The design engineer or project manager is most likely concerned with the veracity and 
flexibility of the program; that is, whether or not the predicted values will accurately 
replicate certain experimental data, and therefore at least reduce the amount of costly 
preliminary testing and experimentation. Exact precision is not always required so long 
as the trend and character of the computer solutions are repetitive and thoroughly doc- 
umented for various missile configurations and flight conditions. The NSWC and 
MISSILE DATCO.M aeroprediction codes have demonstrated a generally comparable 
performance in the prediction of various aerodynamic coefficients. As the prediction 
curves are seen to generally replicate the trend of the experimental data at the correct 
order of magnitude, the NSWC and MISSILE DATCOM codes are found to possess the 
accuracy required to conduct preliminary design and aerodynamic analysis. Introduc- 
tory level research of conventional or arbitrar}’ missile configurations across a wide 
range of input conditions is an ideal application for the .MISSILE DATCO.M prediction 
program. The NSWC aeroprediction code is better suited for a more detailed aerodyna- 
mic investigation of specific missile geometries in support of feasibility studies or inter- 
mediate research. 

An analysis of asymmetric vortex separation is currently not feasible for smaller re- 
search or academic institutions. The study of nonlinear and time variant flow mech- 
anisms has traditionally made use of complex and expensive numerical techniques which 
cannot be operated efficiently by the majority of mainframe computer systems. The re- 
quirement for a low-level, preliminary design code capable of conducting asymmetric 
analysis is evident, and underscores the need to continue development of such programs. 
The VORST.'\B 11 computer program may well represent a first step in fulfilling this 
need. The capabilities and prediction performance of the VORST.AB 11 source code 
should be investigated as this program becomes available to a larger user base. 
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APPENDIX A. BASIC FINNER PITCH DAMPING COEFFICENT 

PREDICTION 




Source: [Ref. 26, p. 113] 
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APPENDIX B. DYNAMIC DERIVATIVE OUTPUT OF PITCH DAMPING 

COEFFICIENT 



THE USAT AUTCUATEO y ISSUE OATOOy • REV 11/BS • 
AEROOYHAMIC umCOS FOR yiSSILE 0QNF1CURATI0N8 
example PROBLEM - PLAMAR RING PUIS TAIL OONF10URATIOM 
STATIC AEWOYNAM1CS FOR BGOY-fiN STT 1 AlO 2 



CASE 1 
PACE 2* 



UACH 

MJIAER 

2.M 



ALTITUJE 

FT 



FLIGHT CXRCITIOMS 

VELOCITY PRESSLJRE TEkTERATLIRE REYTCLDS 

MLAfiEJI 

R/SCC LB/1 R* *2 OCC R 1/H 



SIDESLIP ROLL 
ANGLE ANGLE 

DEO DEC 



REF. 

AREA 

1N««2 

11.R45 



REFERENCE DIMENSIONS 



REF. LENGTH 
LONG. LAT. 

IN IN 

3.7SR 3.7M 



moment ret. center 

L£)NO. VERTICAL 

IN IN 

18.75# #.#•# 



DERIVATIVES (PER DEGREE} 

LONCITUDIHAL LATERAL DIRECTIONAL LONCI TU)IMAL LATERAL DIRECTIONAL 



ALPHA 


CM 


CM 


CA 


CY 


CLN 


CLL 


cna 


CM 


§.#• 


8.888 


8.888 


8.387 


8.808 


8.800 


8 880 


2.568E-81 


-3 424E-81 


4.tf 


1.124 


-1.472 


8 368 


8.808 


8.000 


8.800 


3.834C-#1 


-3.934E-81 


8.## 


2.428 


-3.148 


8.388 


8.808 


8.00# 


8.000 


3.717E-81 


-4.604E-81 


12.6# 


4.184 


-5.327 


8.368 


8.808 


8.800 


8.000 


4.317E-81 


-5.763E-81 


16.8# 


5.882 


-7.78# 


8.369 


8.800 


8 800 


8 800 


4 414C-81 


-6.214E-81 


2# 8# 


7.635 


-18.298 


8.378 


8.000 


8.000 


8 800 


4 194E-81 


-6.17SC-81 


24 88 


9.239 


-12.728 


8.371 


8 800 


8.800 


8 000 


4.83SE-81 


-5.809E-81 


28 8# 


18.883 


-15.826 


8.372 


8.800 


8.00# 


8.800 


4.885C-81 


-5.624E-81 










alpha 


CL 


CO 


CLyCO X-C.P. 












8.8# 


8 800 


8.367 


8.808 -1.323 












4 .8# 


1 .898 


8 445 


2.482 -1.309 












6.8# 


2.353 


8.782 


3.352 -1.297 












12.88 


3.937 


1.213 


3.245 -1.298 












16.88 


5.553 


1 .976 


2.618 -1.323 












28 88 


7.848 


2 959 


2.362 -1.349 












24 88 


6.289 


4.897 


2.823 -1.377 












26.80 


9.416 


5.429 


1.735 -1.363 





CLLB 



Sample of Test Case Header for Flight and Reference Conditions 



THE USAF AUTOMATED MISSILE DATCOM • REV 11/85 . CASE 1 

AEROOYMAMIC METKJOS FOR MISSILE CONFIGURATIONS PACE 3 

EXAi#>LE PR06LEM - PLANAR NINC PLUS TAIL CONFIGURATION 
BODY ♦ 2 fin SETS DYKANIC DERIVATIVES 



uach altitude 

NLA££R 

FT 

2.36 



velocity 
FT/S EC 



- FLIGHT CONDITIONS 

PRESSURE Tti^ERATLRE RTYNOUDS 
KACER 

LB/IN«*2 dec R 1/FT 

3 



SIDESLIP 

ancle 

DEC 

t e# 



ROLL 

ANCLE 

DEG 

• .d« 



REF. 
AREA 
IN**2 
1 1 #43 



REFERENCE DIMENSIONS 



REF, 

LONG. 

IN 

3.75# 



length 
LAT 
IN 

3.75# 



moment ref CENTER 



LONG. 

IN 

16.75# 



vertical 

IN 

e ### 



dynamic derivatives (per DEGREE) 



LPHA 


CNQ 


CHAD 


040KMO 


8 8 


5.564E-81 


3.6e6E-81 


-4.641 19E+08 


4.8 


8.547E-#1 


4.465E-81 


-5 33195C4^ 


a 8 


6 e 22£-#1 


5. 47eE-«t 


-8.5ie44E+O0 


12 8 


9 317E-01 


6.353E-01 


-7.63775E4^88 


16 8 


9.877E-8I 


6.735E-81 


-6.76341E+0# 



Dynamic Derivative Output Including Pitch Damping CoefHcient 
Source; [Ref. 35, p. 123] 
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